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The  performance  effect  of  the  hydrogen  fuel  contaminant  carbon  monoxide  (CO)  for  dry  gas  concen¬ 
trations  ranging  from  1  to  10  ppm  in  proton  exchange  membrane  fuel  cells  (PEMFCs)  is  reported. 
Performance  loss  is  investigated  for  different  operating  conditions  including  several  temperatures,  anode 
relative  humidities,  and  current  densities.  An  analytical  method  is  presented  to  quantify  performance 
loss  due  to  CO  exposure  (A 77),  the  time  required  to  attain  steady  state  of  the  CO  reaction  mechanism 
(brans)*  and  the  effectiveness  of  cell  performance  recovery  in  neat  hydrogen.  For  select  experiments,  the 
magnitude  of  conversion  of  CO  to  carbon  dioxide  (C02)  at  steady  state  (fCo)  is  quantified  with  gas 
chromatography.  The  results  show  (i)  A 7)  is  greatest  at  high  CO  concentrations,  high  current  density,  and 
low  anode  relative  humidity,  (ii)  oxygen  crossover  from  the  cathode  decreases  A  77,  and  (iii)  £co  is  higher 
at  a  temperature  of  60  °C  than  80  °C.  The  latter  may  be  attributed  to  an  increased  amount  of  liquid  water 
in  the  cell  increasing  the  oxygen  crossover.  The  data  also  shows  that  recovery  in  neat  hydrogen  was 
typically  complete.  A  summary  of  all  experimental  results  of  A r\  is  given  in  Appendix  B  as  fitting 
equations  and  parameters  to  assist  in  cell  performance  model  development. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Proton  exchange  membrane  fuel  cells  (PEMFCs)  are  promising 
clean  power  sources  for  transportation  and  stationary  applications. 
When  operated  with  pure  hydrogen  as  the  fuel  and  clean  air  as  the 
oxidant,  PEMFCs  offer  high  efficiencies,  high  power  densities,  and 
zero  emissions  [1,2].  However,  contaminants  in  the  hydrogen  fuel 
stream  (e.g.  CO,  CO2,  H2S,  NH3,  and  organic  compounds)  and/or  in 
the  air  stream  (e.g.  NOx,  SOx,  CO2,  NH3  and  organic  compounds) 
have  been  shown  to  affect  fuel  cell  performance,  sometimes 
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resulting  in  permanent  damage  to  the  membrane  electrode 
assemblies  (MEAs)  [3-5].  Numerous  studies  have  been  reported 
that  characterized  the  effect  of  contaminants  on  cell  performance, 
proposed  reaction  mechanisms  for  contaminants  entering  the  fuel 
cell,  and  described  models  that  predicted  their  effects  on  cell 
performance.  The  majority  of  these  works  focused  on  CO  at 
concentration  levels  between  10  and  100  ppm.  At  these  concen¬ 
trations  significant  performance  loss  was  observed.  However,  based 
on  experiments  conducted  to  support  the  development  of  the  draft 
hydrogen  fuel  quality  specification,  SAE  J2719,  it  is  currently 
believed  that  contaminant  concentration  limits  for  CO  as  low  as 
0.2  ppm  may  be  required  to  maintain  acceptable  cell  performance 
levels  [6]. 

A  significant  detail  in  assessing  the  impact  of  CO  at  concentra¬ 
tions  below  10  ppm  is  the  accurate  determination  of  the  perfor¬ 
mance  loss  in  neat  H2/air  due  to  cell  degradation.  This  loss  is 
expected  to  vary  significantly  with  operating  conditions  but  may 
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also,  as  shown  in  this  work,  vary  between  MEAs  from  the  same  vendor 
and  production  batch.  To  address  this  issue,  we  developed  a  quanti¬ 
tative  methodology  that  separates  the  cell  degradation  from  the 
performance  loss  induced  by  a  contaminant.  The  method  is  applicable 
to  any  operating  condition  with  MEAs  of  any  loading,  batch,  or 
manufacturer.  It  enables  comparison  of  the  contaminant  induced 
performance  impacts  over  a  wide  range  of  operating  conditions. 

In  this  work,  we  report  on  this  method  and  apply  it  to  quantify 
the  performance  impact  of  CO  at  dry  gas  concentrations  of  1,  2,  and 
10  ppm  and  several  operating  conditions  to  affirm  the  proposed 
SAE  J2719  gas  purity  requirements.  The  overpotential  change  of  the 
cell  due  to  CO  exposure  (A tj),  the  time  required  to  attain  steady 
state  of  the  CO  reaction  mechanism  (ttrans)*  and  the  effectiveness  of 
performance  recovery  in  neat  hydrogen  are  reported  and  discussed 
for  each  operating  condition  which  included  various  temperatures, 
relative  humidity,  and  current  densities.  The  role  of  crossover 
oxygen  in  affecting  A rj  and  ttrans  is  also  investigated  and  the 
magnitude  of  CO  to  C02  conversion  when  the  CO  reaction  mecha¬ 
nism  attained  steady  state  is  reported  for  certain  experiments. 


2.  Experimental 

Experiments  discussed  in  this  work  were  performed  at  the 
Hawaii  Fuel  Cell  Test  Facility  (HFCTF),  a  state  of  the  art  laboratory  for 
fuel  cell  research  on  single  cells  and  short  stacks.  HFCTF  features  on¬ 
site  production  of  high-purity  hydrogen,  on-site  air  purification  and 
supply,  and  a  suite  of  diagnostic  equipment.  FCATS™  G050  series  test 
stations  from  Green  Light  Power  Technologies  Inc.  modified  for  on¬ 
line  gas  analysis  at  the  anode  and  cathode  were  used  to  conduct 
experiments  on  50  cm2  single  cell  PEMFCs.  Continuous  1  kHz  resis¬ 
tance  measurements  were  recorded  during  all  polarization  and 
constant  current  experiments  using  a  Model  DOT-1000  Milli- 
Ohmmeter  from  Criterion  Instruments.  A  detailed  description  and 
schematic  of  this  experimental  setup  was  reported  [7].  Careful 
management  of  humidification,  temperature,  and  gas  flow  control 
systems  assured  stable  long-term  operation  of  PEMFCs  [7]. 

Table  1  summarizes  the  fuel  cell  hardware  and  component 
materials  used  in  this  work.  Commercially  available  membrane 
electrode  assemblies  (MEA)  with  0.375  mg  Pt  cm-2  loading  and  50% 
Pt/C  catalysts  were  purchased  from  Ion  Power.  Membrane  and 
electrode  thicknesses  were,  40  pm  and  20  pm,  respectively,  and 
were  determined  from  MEA  cross  sections  via  scanning  electron 
microscopy  (SEM).  The  cells  were  assembled  using  25  BC  gas 
diffusion  layers  (GDL)  from  SGL  Technologies  on  both  electrodes. 
Standard  protocols  were  employed  for  cell  assembly,  conditioning, 
and  diagnostics.  Prior  and  subsequent  to  contaminant  testing, 
beginning  of  test  (BOT)  and  end  of  test  (EOT)  diagnostics  were 
performed  on  the  cell.  Diagnostics  included  polarization  curves 
with  H2/air  and  H2/02  (VI  curves). 

Table  2  summarizes  the  operating  conditions  used  in  this  work. 
H2/air  and  H2/02  polarization  curves  were  obtained  using  ultra- 

Table  1 

Cell  hardware  and  components. 


high  purity  gases  humidified  prior  to  flowing  into  the  anode  and 
cathode  of  the  cell.  Polarization  curves  were  measured  in  current 
control  mode,  starting  at  a  current  density  of  1.2  A  cm-2  with 
0.2  A  cm-2  steps  towards  open  circuit  (OC).  Cells  were  operated  at 
a  flow  stoichiometry  of  two  at  the  anode  and  cathode  for  each 
current  density  when  air  was  the  oxidant  gas.  For  curves  measured 
with  02,  the  flow  rates  were  the  same  as  those  in  the  air  curves 
resulting  in  a  cathode  stoichiometry  of  9.5.  An  additional  data  point 
was  recorded  at  0.1  A  cm-2.  Anode  and  cathode  flow  rates  at  this 
measurement  point  were  equivalent  to  the  flows  at  a  current 
density  of  0.2  A  cm-2.  At  each  point,  the  current  was  held  constant 
for  15  min.  The  voltage  values  reported  were  averaged  over  the  last 
5  min.  The  cell  was  subsequently  held  for  1  min  at  OC  and  the 
maximum  recorded  voltage  was  used  as  the  open  circuit  voltage. 

Long  duration  constant  current  density  experiments  were  con¬ 
ducted  at  fixed  operating  conditions  to  enable  the  quantification  of 
the  magnitude  of  performance  loss  due  to  CO  exposure.  The  voltage 
change  of  cells  operating  continuously  in  current  control  mode  for 
durations  of  100-200  h  with  neat  H2/air  was  measured  at  different 
conditions  to  characterize  the  performance  loss  of  the  MEA. 

For  CO  exposure  experiments,  the  cell  was  operated  for  a  period 
of  time  with  neat  hydrogen  and  air  followed  by  a  period  where 
a  dry  gas  mixture  of  hydrogen  and  CO  was  continuously  injected 
into  the  humidified  neat  hydrogen  stream.  The  specified  humidi¬ 
fication  was  maintained  by  increasing  the  temperature  of  the 
humidifier  to  compensate  for  the  dry  contaminant  injection. 

3.  Methodology 

This  section  describes  how  fco*  A)],  and  ttrans  were  calculated. 
Values  for  ?Co  were  quantified  from  CO  and  C02  concentrations 
measured  using  gas  chromatography  at  the  anode  inlet  and  outlet. 
The  molar  flow  rates  of  CO  and  C02  entering  and  exiting  the  fuel  cell 
were  quantified  following  the  equations  given  in  [7]  and  used  to 
determine  the  quality  of  closure  for  the  carbon  balance  (QcBai)-  QcBai 
was  calculated  using  Eq.  (1),  where  hi  out  and  hl  in  represent  the 
total  molar  flow  rate  of  carbon  containing  impurity  species  out  of 
and  into  the  electrode,  respectively.  A  value  of  one  indicates 
complete  closure  of  the  molar  flow  balance  of  carbon  at  the  anode. 


QcBai 


H-i,  out 
_i _ 

Y2  fy, in 
i 


(1) 


The  conversion  of  CO  to  C02  at  steady  state  was  calculated  using 
Eq.  (2).  The  calculation  assumes  closure  of  the  mass  balance,  which 
was  validated  in  previous  work  where  the  C  mole  balance  closed  to 
within  3%  [7].  fco  was  calculated  as  shown  in  Eq.  (2),  where  fco  was 
the  conversion  of  CO,  hCOin  and  hCoout  were  the  molar  flow  rates  of 
CO  entering  and  exiting  the  cell,  respectively. 


?co 


t*COin  ^COout 

r*COin 


(2) 


Hardware 

General 

Anode 

Cathode 

MEA 

Ion  power 

Area 

50  cm2 

Membrane  thickness 

40  pm 

Electrode  thickness 

20  pm 

20  pm 

Flow-field 

2  Channel  serpentine 

3  Channel  serpentine 

Catalyst 

50%  Pt/C 

50%  Pt/C 

Loading 

0.375  mg  Ptcm-2 

0.375  mg  Ptcm-2 

GDL 

SGL-BC25 

SGL-BC25 

Gaskets 

178  pm,  Teflon 

203  pm,  Teflon 

Table  2 

Operating  conditions. 


Operating  conditions 

General 

Anode 

Cathode 

Relative  humidity 

100/80/60% 

50% 

Stoichiometry 

2 

2 

Fuel/oxidant 

h2 

Air/oxygen/hydrogen 

Pressure 

Temperature 

Current  Density 

48.3  l<Pag 

80/60  °C 

1/0.6/0.2  A  cm-2 

48.3  kPag 
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This  section  describes  the  experimental  and  analytical  meth¬ 
odology  developed  to  quantify  the  impact  of  trace  level  contami¬ 
nants  on  fuel  cell  performance.  The  methodology  separates  the 
performance  loss  that  occurs  during  operation  with  neat  H2/air 
from  that  caused  by  the  contaminant.  It  also  provides  specific 
procedures  for  determining  At?  and  ttrans- 

Fig.  1  shows  the  measured  cell  voltage  data  (black,  left  axis)  and 
At?  (grey,  right  axis)  for  an  experiment  conducted  at  constant 
current  with  1  ppm  CO.  Conditions  are  given  within  the  figure.  The 
experiment  is  divided  into  three  distinct  phases  which  are  labeled 
in  Fig.  1. 

•  Phase  1  is  the  pre-poisoning  period  conducted  with  neat  H2/air. 
The  performance  decay  during  this  phase  is  used  to  extrapolate 
the  projected  performance  decay  of  the  cell  under  neat  H2/air 
operation  for  a  time  equivalent  to  the  duration  of  the  entire 
experiment  (i.e.  Phases  1  through  3).  The  length  of  Phase  1 
required  to  achieve  accurate  extrapolation  is  approximately  1  /3 
of  the  length  of  the  entire  experiment.  The  quantification  of  the 
fitting  interval  is  discussed  in  Appendix  A.  Performance  loss 
during  Phase  1  is  attributed  to  changes  in  the  MEA. 

•  Phase  2  begins  with  the  introduction  of  contaminant  into  the 
anode  feed  stream.  In  the  experiments  reported  in  this  work, 
contaminant  concentration  was  held  constant  throughout 
Phase  2.  Performance  loss  during  this  phase  includes  the  losses 
due  to  MEA  effects  which  were  estimated  by  extrapolation  of 
the  Phase  1  data  plus  the  loss  due  to  the  contaminant.  Perfor¬ 
mance  loss  occurring  during  operation  with  pure  gases  was 
assumed  to  be  continuous  and  independent  of  losses  resulting 
during  exposure  to  contaminants. 

•  Phase  3,  the  recovery  period,  begins  when  the  flow  of  the 
contaminant  is  stopped  and  operation  on  neat  H2/air  is 
resumed.  Performance  changes  during  this  phase  include 
continuing  losses  due  to  degradation  of  the  MEA  plus  the 
performance  recovery  resulting  from  the  reaction  of  CO  to  C02 
and  CO  desorption. 

Analysis  of  the  data  from  all  three  phases  allows  the  perfor¬ 
mance  impact  due  to  contaminant  exposure  to  be  accurately 
quantified  by  separating  it  from  the  performance  decrease  due  to 
MEA  degradation.  It  also  enables  accurate  determination  of  the 


Time  [hours] 

Fig.  1.  Top:  experimentally  measured  cell  data  versus  time  for  exposure  to  dry  gas 
concentration  of  1  ppm  CO  at  45  h.  Extrapolated  fit  of  pre-poisoning  period  throughout 

experimental  duration  shown  as  dashed  line  ( - ).  Bottom:  calculated  values  for 

overpotential  change  (At?)  versus  time  due  to  CO  exposure. 


extent  of  recovery  following  removal  of  the  contaminant  from  the 
fuel  stream.  All  contaminant  experiments  were  analyzed  in  this 
manner. 

Fig.  1  illustrates  the  application  of  the  methodology  described 
above.  The  top  line  shows  the  measured  cell  voltage  response  while 
the  bottom  line  shows  the  calculated  change  in  cell  overpotential 
attributed  to  the  CO,  i.e.  At?.  Phase  1,  lasted  from  0  to  45  h  which 
represented  approximately  1  /3  of  the  duration  of  the  experiment 
and  showed  the  measured  cell  voltage  under  neat  H2/air  operation. 
The  Phase  1  data  was  fitted  using  Eq.  (3)  to  extrapolate  the 
performance  degradation  through  the  test  time  of  145  h.  Eq.  (3)  is 
purely  empirical  and  consists  of  an  exponential  and  linear 
component  since  it  was  determined  to  achieve  the  best  overall  fit  to 
the  data  from  several  different  types  of  fitting  equations.  The 
extrapolated  data  is  represented  by  the  dashed  line  ( — )  in  the  top 
graph  in  Fig.  1. 

Vfit(t)  =  Pi  'exp(“^“)  +P3  +  t  p4  (3) 

with  pi  >  0,  p2  >  0,  P3  >  0,  P4  <  0 

Although  the  projected  fit  and  the  cell  performance  after 
recovery  are  in  good  agreement,  it  is  important  to  note  that  this 
‘front-loaded’  fitting  procedure  relied  entirely  on  data  prior  to 
contaminant  exposure.  Results  at  our  facility  and  other  research 
institutions  have  shown  that  performance  recovery  in  neat  H2/air 
subsequent  to  exposure  to  a  contaminant  may  remain  incomplete 
[8-10].  Thus,  performance  at  the  end  of  the  recovery  phase  (Phase 
3)  cannot  be  used  for  data  fitting.  For  each  experiment,  the  reli¬ 
ability  of  the  ‘front-loaded’  fit  was  qualitatively  assessed  using 
three  criteria:  (i)  there  was  no  significant  deviation  of  At?  from  zero 
over  the  pre-poisoning  period,  (ii)  the  calculated  overpotentials  did 
not  contain  any  negative  values,  which  would  indicate  an  over¬ 
compensation  of  the  performance  degradation  of  the  fuel  cell  by 
the  fit,  and  (iii)  during  the  final  stages  of  cell  recovery  a  constant 
overpotential  was  observed. 

During  phase  2,  from  45  to  100  h,  1  ppm  of  CO  was  continuously 
supplied  to  the  fuel  inlet  stream.  Upon  introduction  of  the  CO,  the 
rate  of  loss  in  cell  performance  increased  significantly  due  to  CO 
adsorption  reducing  the  active  Pt  reaction  sites  available  for  the 
hydrogen  oxidation  reaction  (FIOR).  As  the  reactions  with  CO  that 
were  involved  in  catalyst  poisoning  attained  steady  state,  the 
measured  rate  of  performance  loss  became  identical  to  the  rate  of 
performance  loss  from  the  extrapolated  fit  of  Phase  1.  The 
measured  cell  voltage  data  in  Fig.  1  indicated  that  steady  state 
conditions  were  reached  at  approximately  70  h  and  were  main¬ 
tained  until  the  contaminant  flow  was  stopped,  at  100  h  of 
operation. 

Phase  3,  the  recovery  phase  was  initiated  at  100  h  by  stopping  the 
flow  of  the  contaminant  and  returning  to  operation  with  neat  H2/air 
which  was  continued  until  the  experiment  was  terminated  at  145  h. 
During  this  phase,  the  cell  performance  initially  recovered  rapidly  and 
then  asymptotically  approached  the  MEA  performance  predicted  by 
the  extrapolated  fit  of  Phase  1.  This  indicated  that  the  performance 
recovery  at  these  operating  conditions  was  nearly  complete. 

Overpotential  change  (At?)  data  was  calculated  by  subtracting  the 
measured  cell  voltage  from  the  extrapolated  neat  H2/air  perfor¬ 
mance  (Eq.  4).  The  At?  data  is  represented  by  the  lower  line  in  Fig.  1. 

Mti)  =  Vfit (tf)  -  Vc ell(tf)  (4) 

During  Phase  1,  under  neat  H2/air  operation,  the  calculated  At? 
was  zero,  as  expected.  Following  introduction  of  the  CO  into  the 
hydrogen  stream  in  Phase  2,  the  At?  increased  rapidly,  and 
approached  a  steady  state  value  of  25  mV.  It  is  important  to  note 
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Fig.  2.  Example  of  quantitative  analysis  of  processed  contaminant  data  to  determine 
ttrans.  the  time  period  until  steady  state  conditions  are  reached  and  A?/,  the  over¬ 
potential  change  due  to  contamination. 


that  the  performance  degradation  due  to  CO  exposure  is  similar  in 
magnitude  to  the  projected  neat  hydrogen  loss  over  the  same 
period  of  time.  This  emphasized  the  need  for  a  methodology  that 
separates  the  performance  impact  of  CO  from  the  performance 
degradation  due  to  neat  FU/air  operation. 

We  also  developed  a  procedure  to  estimate  the  time  required  for 
the  poisoning  process  to  reach  steady  state,  which  we  refer  to  as  the 
transient  time  of  the  poisoning  process  ttrans-  To  facilitate  compar¬ 
ison  of  the  overpotential  data  from  experiments  with  different 
times  before  poisoning  (Phase  1)  overpotential  data  are  shown  so 
that  the  contaminant  injection  appears  to  begin  at  10  h.  Fig.  2 
shows  A?y  for  the  experiment  shown  in  Fig.  1  plotted  in  this 
manner.  A  five-step  process  was  developed  to  reproducibly  esti¬ 
mate  ttrans  and  Ary.  The  methodology  is  illustrated  in  Fig.  2  using  the 
data  of  the  contaminant  experiment  shown  in  Fig.  1 : 

(1)  Identify  an  arbitrary  time,  tarb,  that  is  near  the  beginning  but 
clearly  within  the  steady  state  region.  In  this  region,  Ary  typi¬ 
cally  exhibits  a  linear  progression  with  a  slope  near  zero. 


(2)  Apply  a  linear  fit  from  tarb  to  the  end  of  the  contaminant 
injection,  tstop. 

(3)  Plot  the  fit  and  the  change  of  overpotential  data  into  the  same 
graph  and  identify  the  time,  tss,  when  the  value  of  Ary  and  the 
linear  fit  begin  to  diverge.  This  time  indicates  when  steady 
state  has  been  reached. 

(4)  Subtract  the  time  of  the  contaminant  injection,  qn,  from  tss  to 
determine  the  transient  time,  ttrans.  required  to  reach  steady 
state. 

(5)  Average  the  Ary  data  between  ttrans  and  tstop  to  determine  Ary  at 
steady  state. 

In  the  data  shown  in  Fig.  2,  the  rate  of  change  in  Ary  significantly 
slowed  and  approached  a  slope  of  0  at  approximately  30  h.  This 
indicated  the  approximate  onset  of  steady  state.  The  arbitrary  time 
tarb  was  selected  to  be  35  h  to  assure  the  fit  would  be  well  within  the 
steady  state  region.  The  data  from  tarb  (35  h)  to  the  end  of  the 
contaminant  injection  (tstop  =  65  h)  was  fitted  using  a  linear 
regression.  The  fit  was  represented  by  the  dashed  horizontal  line 
shown  in  Fig.  2.  Deviation  of  the  fit  from  the  Ary  data  indicated  that 
steady  state  conditions  were  reached  at  31  h  and  thus  ttrans  was  21  h. 
The  change  in  overpotential  due  to  CO  exposure  at  steady  state  was 
determined  by  averaging  Ary  from  31  to  65  h  and  found  to  be  25  mV. 
Values  for  Ary  for  CO  concentrations  between  1  and  10  ppm  at  several 
conditions  are  presented  in  the  following  section.  These  results  are 
expected  to  be  of  significant  value  for  the  determination  of  reaction 
rate  constants  and  contaminant  model  verification. 

4.  Results  and  discussion 

In  this  section,  we  present  data  and  discuss  the  effects  of  CO 
concentrations  of  1, 2,  and  10  ppm  at  different  operating  conditions 
on  Ary,  ttrans,  ?Co.  and  the  performance  recovery  in  neat  H2/air. 
Unless  otherwise  noted  CO  concentrations  are  reported  as  dry  gas 
concentrations.  Experiments  discussed  in  this  section  were  per¬ 
formed  with  PEMFCs  operated  at  the  conditions  given  in  Table  3, 
which  also  provides  a  summary  of  all  quantitative  results  including 
Ary,  ttrans.  QcBai.  and  ?co-  Data  were  analyzed  using  the  methodology 
described  in  Section  3. 

4.J.  Effect  of  CO  on  cell  overpotential  and  transient  time 

Fig.  3  shows  Ary  versus  time  for  cells  exposed  to  2  ppm  CO  at 
current  densities  of  0.2, 0.6,  and  1.0  A  cm-2,  at  a  cell  temperature  of 


Table  3 

Transient  time  ttrans  and  overpotential  change  Ary  for  1,  2,  and  10  ppm  CO  at  various  operating  conditions. 


[CO]  [ppm] 

ECA  [m2  gpt1] 

Temp.  (°C)/Anode  R.H.  (%) 

Current  density  [A  cm  2] 

Cathode  gas 

hxans[h] 

At/  [mV] 

0cBal[%] 

£co[%] 

2 

84.07 

60/100% 

1 

Air 

14 

246 

98.1 

70.3 

2 

71.87 

60/100 

0.6 

Air 

29 

159 

- 

- 

2 

85.12 

60/100 

0.2 

Air 

69 

19 

- 

- 

2 

72.32 

60/100 

1 

o2 

n/a 

35 

- 

- 

2 

73.28 

60/100 

1 

h2 

n/a 

300 

- 

- 

1 

65.12 

80/100 

1 

Air 

25 

25 

- 

- 

2 

69.28 

80/100 

1 

Air 

14 

45 

- 

- 

10 

69.33 

80/100 

1 

Air 

9 

241 

110.9 

33.4 

1 

88.89 

80/80 

1 

Air 

22 

25 

- 

- 

2 

80.25 

80/80 

1 

Air 

13 

48 

- 

- 

10 

83.88 

80/80 

1 

Air 

13 

238 

106.5 

21.1 

1 

85.94 

80/60 

1 

Air 

36 

27 

- 

- 

2 

71.98 

80/60 

1 

Air 

13 

77 

96.9 

43.7 

10 

83.47 

80/60 

1 

Air 

12 

323 

102.3 

32.8 

1 

60.94 

60/100 

1 

Air 

29 

173 

100.6 

61.7 

10 

85.84 

60/100 

1 

Air 

8 

343 

97.7 

91.3 
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60  °C,  and  an  anode  relative  humidity  of  100%.  It  should  be  noted 
that  the  recovery  phase  in  pure  H2  was  not  performed  for  the 
0.2  A  cm-2  experiment.  The  results  showed  a  significant  decrease  of 
brans  with  increasing  current  density.  Compared  to  a  current  density 
of  1.0  A  cm-2,  ttrans  was  2.1  and  4.9  times  greater  for  the  0.6  and 
0.2  A  cm-2  experiments,  respectively.  Because  the  flow  stoichiom¬ 
etry  and  injection  concentration  of  CO  were  held  constant,  the 
absolute  flow  rate  and  thus  the  contaminant  molar  flow  rate  into  the 
cell  was  reduced  by  a  factor  of  1.7  and  5.0,  for  the  0.6  and  0.2  A  cm-2 
experiments,  respectively.  This  suggested  a  directly  proportional 
relationship  between  the  molar  flow  rate  of  CO  into  the  cell  and  ttrans- 

A  significantly  greater  performance  loss  was  observed  at  higher 
current  densities.  The  value  for  A rj  at  the  current  densities  of  0.2, 
0.6,  and  lAcm-2  was  19,  159,  and  246  mV,  respectively.  This 
increase  of  A  rj  with  current  density  likely  resulted  from  the 
competitive  adsorption  of  CO  and  hydrogen.  At  higher  current 
densities  a  greater  fraction  of  the  electrochemically  active  surface 
area  was  required  to  sustain  the  hydrogen  oxidation  reaction, 
which  in  the  presence  of  CO,  decreases  cell  voltage.  For  these 
experiments,  the  exact  Pt-CO  or  Pt-H  surface  coverage  was  not 
known;  however,  the  trend  of  increased  overpotential  for  a  fixed  CO 
concentration  at  higher  operating  current  densities  was  consistent 
with  the  modeling  predictions  of  previous  reports  for  higher  CO 
concentrations  (>10  ppm)  [3,4,11]. 

4.2.  The  role  of  crossover  oxygen  on  the  oxidation  of  CO 

Some  of  the  CO  that  adsorbs  to  the  Pt/C  catalyst  within  the  fuel 
cell  reacts  to  form  C02.  Various  studies  have  reported  that  the 
reaction  of  CO  to  form  C02  may  proceed  via  a  chemical  reaction,  an 
electrochemical  reaction,  or  a  combination  of  the  two.  A  generally 
accepted  mechanism  for  the  electrochemical  oxidation  of  CO  that 
was  linearly  adsorbed  onto  Pt  is  given  by  Eqs.  (5)  and  (6)  [4],  [11  — 
13]. 

CO  +  Pt  Pt-COads  (5) 

Pt  -  COads  +  H20— >Pt  +  C02  +  2H+  +  2e~  (6) 

A  similar  mechanism  has  been  reported  for  the  electrochemical 
oxidation  of  CO  that  was  bridge  bonded  on  Pt  [12].  In  either  of  these 


Fig.  3.  Comparison  of  overpotentiai  change  (At/)  versus  time  for  ceils  exposed  to  dry 
gas  concentration  of  2  ppm  CO  at  60  °C;  0.2,  0.6,  and  1  A  cm-2. 


cases  water  was  considered  to  be  the  source  of  oxygen  for  the 
electrochemical  oxidation  of  CO. 

Literature  has  indicated  that  the  electro-oxidation  of  CO 
requires  a  potential  of  at  least  350  mV  on  a  platinum/ruthenium 
alloy  catalyst  versus  a  reversible  hydrogen  electrode  (RHE)  [14-16] 
to  occur.  For  fuel  cells  using  Pt/C  as  the  catalyst,  such  as  those  in  this 
work,  the  ignition  potential  for  CO  electro-oxidation  is  expected  to 
be  even  higher.  As  shown  in  Table  3,  the  maximum  value  for  A  rj 
measured  in  this  work  was  343  mV  for  a  10  ppm  CO  exposure  at 
60  °C  and  an  anode  relative  humidity  of  100%.  Most  of  the  A  rj  values 
calculated  were  substantially  lower.  Previous  work  has  also  shown 
that  the  potential  for  the  onset  of  the  CO  electro-oxidation  and  its 
peak  potential  are  dependent  on  cell  temperature  and  cell 
humidification.  Both  were  reported  to  decrease  with  increasing  cell 
temperature  [14,15].  Furthermore,  the  CO  stripping  peak  potential 
was  shown  to  increase  when  the  working  electrode  of  the  cell  was 
operated  with  a  less  than  fully  humidified  gas  feed  stream  [16]. 
Therefore,  the  overpotentials  calculated  in  this  work  are  likely 
below  the  potential  onset  of  CO  electro-oxidation,  and  are  thus 
unlikely  to  occur. 

A  second  possible  mechanism  for  CO  oxidation  is  a  chemical 
reaction  with  oxygen  which  enters  the  anode  by  diffusing  across 
the  MEA  from  the  cathode.  The  impact  of  crossover  oxygen  on  the 
oxidation  of  platinum  adsorbed  CO  has  gained  importance  in  recent 
years  since  the  thickness  of  proton  exchange  membranes  in  fuel 
cells  has  significantly  decreased  from  180  to  25  pm  which  has 
thereby  increased  the  rate  of  oxygen  crossover.  Evidence  of  oxygen 
crossover  and  reaction  with  CO  was  given  by  Zhang  et  al.  [17],  who 
reported  that  reaction  of  CO  to  C02  may  occur  as  a  catalytic  surface 
redox  mechanism  such  as  that  shown  by  Eqs.  (7)  and  (8),  or  through 
the  formation  of  adsorbed  hydroxyls,  OFIads,  that  react  further  with 
surface  bound  CO  electrochemically  as  shown  in  Eqs.  (9)  and  (10). 


02+2Pt~2(Pt-0ads)  (7) 

Pt  —  Oads  +  Pt  —  COads  — ►  2Pt  +  C02  (8) 

Pt  -  Oads  +  Pt  -  Hads  <-  Pt  -  OHads  (9) 

Pt  -  OHads  +  Pt  -  COads— >2Pt  +  C02  +  H+  +  e~  (10) 


To  further  investigate  the  role  of  oxygen  as  a  possible  reactant 
with  CO  at  low  concentrations,  experiments  were  performed  with 
different  cathode  gas  compositions  which  included  100%  hydrogen 
or  100%  oxygen.  Fig.  4  shows  the  change  in  cell  overpotentiai  versus 
time  for  a  cell  operated  at  60 °C,  lAcm-2,  an  anode  relative 
humidity  of  100%  and  2  ppm  CO  in  the  fuel  stream.  The  cathode 
inlet  humidity  was  100%  when  hydrogen  was  supplied  to  the 
cathode  and  the  cell  was  operated  in  an  electrically  driven  mode. 
The  cathode  inlet  humidity  was  50%  when  air  or  02  was  supplied  to 
the  cathode  and  the  cells  were  operated  galvanically.  The  values  for 
At/  were  300, 250,  and  35  mV,  when  the  cathode  gas  was  hydrogen, 
air,  or  oxygen,  respectively.  Clearly,  the  presence  and  concentration 
of  oxygen  at  the  cathode  significantly  reduced  At/.  Oxygen  crossing 
over  from  the  cathode  to  the  anode  presumably  reacted  with  the 
CO  to  form  C02.  Based  on  the  At/  values  we  believe  that  the  reaction 
was  a  chemical  process  such  as  that  given  by  Eqs.  (7)  and  (8).  More 
work  is  required  to  characterize  and  quantify  the  impact  of  cross¬ 
over  oxygen  to  gain  a  greater  understanding  of  the  how  oxygen  on 
the  cathode  affects  the  performance  of  a  PEMFC  with  CO  in  the 
anode  feed,  specifically  at  low  CO  contaminant  concentrations. 

Fig.  5  shows  At/  data  calculated  at  steady  state  with  the  method 
discussed  above  for  1,  2,  and  10  ppm  dry  CO  concentration  expo¬ 
sures  at  several  operating  conditions.  To  allow  for  more  accurate 
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Fig.  4.  Comparison  of  overpotential  change  (At/)  for  cells  exposed  to  dry  gas  concen¬ 
tration  of  2  ppm  CO  at  60  °C,  and  (i)  H2/H2)  (ii)  H2/air,  and  (iii)  H2/02  anode/cathode 
reactant  gases. 

comparisons  the  contaminant  concentrations  were  calculated  to 
account  for  the  amount  of  water  vapor  in  the  inlet  gas  stream. 
Operating  conditions  included  a  60  °C  cell  temperature  with  an 
anode  inlet  humidity  of  100%,  and  a  cell  temperature  of  80  °C  with 
anode  inlet  humidities  of  100, 80,  and  60%.  The  results  showed  that 
At/  increased  with  CO  concentration  at  all  operating  conditions. 
This  confirmed  that,  as  previously  discussed,  higher  concentrations 
of  CO  generally  resulted  in  greater  performance  losses.  The  values 
for  Ary  were  lower  at  80  °C  than  at  60  °C.  This  may  have  indicated 
that  the  CO  coverage  on  the  Pt  catalyst  is  higher  at  lower  operating 
temperatures  due  to  thermal  kinetic  effects.  However,  the  partial 
pressure  of  water  was  27.4  kPa  lower  at  60  °C  which  increased  the 
actual  concentration  of  CO  in  the  feed  stream  and  may  have 
increased  the  equilibrium  surface  coverage.  Water  transport  within 
the  cell  may  have  affected  the  transport  of  oxygen,  which  as  shown 
in  Fig.  4,  mitigates  performance  loss. 


Fig.  5.  Comparison  of  overpotential  change  (At/)  at  steady  state  for  cells  exposed  to  dry 
gas  concentrations  of  1,  2, 10  ppm  CO  at  several  operating  conditions. 


Fig.  6  was  included  to  emphasize  how  the  mass  transport 
limitations  changed  with  operating  conditions  for  the  particular 
cell  design,  GDL,  and  MEA  used  in  this  study.  The  figure  shows  fuel 
cell  performance  curves  and  1  kHz  resistances  measured  at  80  and 
60  °C  with  100%  humidity  at  the  anode  inlet  and  50%  humidity  at 
the  cathode  inlet.  No  change  in  the  cell  resistance  was  observed  at 
80  or  60  °C;  however,  a  significant  performance  deviation  was 
apparent  at  current  densities  above  40  mA  cm-2.  This  indicated 
that  the  performance  at  60  °C  was  impacted  by  higher  cathode 
mass  transport  limitations.  Considering  Figs.  5  and  6  together,  the 
results  suggested  that  the  physical  state  of  water  within  the  cell 
may  have  also  affected  the  cell  performance  when  CO  was  present. 

Fig.  5  also  shows  A rj  values  of  cells  that  were  operated  with 
different  anode  humidities.  The  data  did  not  show  a  clear  trend 
with  humidification.  At  CO  concentrations  near  1  ppm  A  rj  was 
similar  for  all  anode  humidification  conditions.  At  higher  CO 
concentrations  a  humidity  of  60%  resulted  in  the  largest  perfor¬ 
mance  loss  and  the  lowest  performance  loss  was  observed  at  80%. 
The  exact  reasons  for  the  different  values  of  At/  are  not  known  since 
many  of  the  processes  that  take  place  within  the  fuel  cell  are 
interlinked  which  makes  it  difficult  to  clearly  identify  trends  of  At/ 
with  certain  operating  parameters.  In  general  though,  the  trends 
observed  here  with  At/  and  operating  conditions  are  consistent 
with  those  reported  previously  for  contaminant  concentrations 
>10  ppm  CO  [3,4,7,18]. 

Fig.  7  shows  results  for  the  conversion  of  CO  to  CO2,  ?co.  deter¬ 
mined  by  gas  chromatography  at  steady  state  conditions  for 
selected  experiments.  The  CO  concentration  values  in  the  figure 
account  for  the  water  pressure  in  the  humidified  gas  stream.  For 
60  °C,  ?co  increased  with  CO  concentration.  This  trend  was  similar 
to  that  observed  in  Fig.  5  for  At/.  It  indicated  that  higher  concen¬ 
trations  were  likely  related  to  an  increased  rate  of  reaction  due  to 
higher  Pt-CO  catalyst  coverage.  The  conversion  of  CO  was  consis¬ 
tently  lower  at  80  °C.  This  may  be  directly  related  to  a  lower  catalyst 
coverage  of  CO  at  the  higher  temperature  and/or  indirectly  related 
to  higher  amounts  of  liquid  water  at  60  °C.  According  to  Basura 
et  al.  the  presence  of  liquid  water  enhances  the  transport  of  oxygen 
across  the  MEA  [19],  and  according  to  Zhang  et  al.,  the  presence  of 
oxygen  at  the  anode  increases  the  chemical  oxidation  of  CO  to  CO2 
[17].  Both  together  would  consequently  result  in  an  increased  fco- 

Other  literature  shows  that  the  oxygen  crossover  rate  is  also 
dependent  on  the  operating  conditions  of  the  fuel  cell.  The  oxygen 


An/Cath:  48.3/48.3  kPag,  2/2stoich,  100/50%  RH,  HJAir 


Fig.  6.  VI  curves  and  1  kHz  resistances  for  fuel  cells  operated  at  80  and  60  °C  with  an 
anode/cathode  humidity  inlet  humidity  of  100/50%. 
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Fig.  7.  Magnitude  of  CO  to  C02  conversion  quantified  from  gas  chromatography  data  at 
different  operating  conditions. 


crossover  rate  was  reported  to  increase  with  temperature,  with  the 
relative  humidity  of  the  fuel  cell,  and  with  the  cathode  oxygen 
concentration  [19-22].  Ma  et  al.,  suggested  that  the  permeability  of 
oxygen  in  Nation  is  mainly  influenced  by  thermal  diffusion  [20], 
while,  as  mentioned  above,  Basura  et.al  concluded  that  the  oxygen 
diffusion  coefficient  is  related  to  the  water  content  of  the 
membranes  and  suggested  that  the  aqueous  phase  is  dominantly 
involved  in  the  diffusion  pathway  [19].  A  comparison  of  the  results 
in  Figs.  6  and  7  may  confirm  the  literature  results  since  they 
respectively  show  an  impact  of  liquid  water  on  the  cell  performance 
and  an  increase  of  the  conversion  of  CO  to  CO2. 


4.3.  Self-induced  recovery 

Fig.  8  shows  data  for  the  transition  from  phase  2,  the  CO  exposure 
phase,  to  phase  3,  the  recovery  phase  in  neat  H2  and  air.  The  time 


Fig.  8.  Comparison  of  overpotentiai  change  (A rj)  versus  recovery  time  for  dry  CO  gas 
concentrations  of  1,  2,  and  10  ppm,  60  °C  and  100%  relative  humidity. 
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Fig.  9.  Comparison  of  overpotentiai  change  (At/)  versus  recovery  time  for  dry  CO  gas 
concentration  of  2  ppm  and  several  operating  conditions  which  are  given  within  the 
figure. 
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scale  in  the  figure  is  shifted  so  that  the  recovery  phase  starts  at  0  h. 
The  data  indicated  a  decrease  in  A r/  versus  time  when  switching 
from  1,  2,  and  10  ppm  CO  exposure  to  neat  H2/air  operation. 
Although  A ri  at  steady  state  varied  from  173  to  343  mV  as  CO 
concentration  increased  from  1  to  10  ppm,  the  recovery  of  the  three 
cell  voltages  is  nearly  identical.  During  the  first  hour  of  recovery  At/ 
decreased  significantly  and  then  asymptotically  approached  values 
that  were  only  a  few  mV  above  zero  over  several  tens  of  hours.  This 
indicated  nearly  complete  performance  recovery  for  these  concen¬ 
trations  and  operating  conditions  without  an  externally  driven 
recovery  process.  The  shape  and  similarity  of  the  A?/  plots  was 
indicative  of  rapid  CO  desorption  from  the  catalyst  surface  upon 
reintroduction  of  the  neat  hydrogen.  The  time  scale  for  performance 
recovery  in  pure  H2  appeared  to  be  independent  of  the  CO  concen¬ 
tration  at  these  operating  conditions  and  concentrations. 

Fig.  9  shows  this  self-induced  recovery  for  a  variety  of  operating 
conditions  following  poisoning  with  2  ppm  CO.  Conditions  include 
cell  current  densities  of  0.6  and  1.0  A  cm”2  at  60  °C,  and  anode/ 
cathode  relative  humidity  of  60/50,  80/50,  and  100/50%  at  80  °C. 
Again  the  time  scale  was  shifted  so  that  the  recovery  phase  started 
at  0  h.  The  data  indicated  two  trends  with  operating  conditions: 
lower  cell  temperature  and  lower  anode  humidification  resulted  in 
a  slower  recovery  process.  The  temperature  effect  was  expected, 
since  the  desorption  rate  is  strongly  temperature  dependent  and 
increases  with  increasing  temperature.  At  lower  humidification, 
however,  while  considering  similar  conversion  of  CO,  as  shown  in 
Fig.  9,  a  reduced  activity  of  water  may  have  contributed  to  this 
effect.  In  any  case,  the  data  indicated  that  the  operating  conditions 
tested  in  this  work  did  not  have  a  significant  effect  on  the 
completeness  of  performance  recovery.  All  cells  attained  significant 
performance  recovery  within  the  first  five  hours  of  being  subjected 
to  neat  hydrogen. 


5.  Conclusions 

An  analytical  methodology  was  developed  that  enabled  accurate 
quantitative  determination  of  the  performance  loss  associated  with 
the  presence  of  low-level  CO  contaminants  (<10  ppm)  within  the 
anode  feed  stream  of  a  PEMFC.  The  methodology  enabled 
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comparison  of  the  contaminant  impact  between  different  MEAs, 
MEAs  operated  at  various  conditions,  and  MEAs  exposed  to  different 
contaminant  concentrations.  The  methodology  divided  the  experi¬ 
ment  into  three  distinct  phases:  a  pre-poisoning  phase,  a  poisoning 
phase,  and  a  recovery  phase  with  the  minimum  duration  of  the  pre¬ 
poisoning  being  approximately  1/3  of  the  entire  experiment  dura¬ 
tion.  Extrapolation  of  the  pre-poisoning  data  then  allowed  accurate 
prediction  of  the  neat  H2/air  performance  degradation  of  the  cell.  It 
also  enabled  us  to  separate  and  accurately  determine  the  degrada¬ 
tion  due  to  contaminant  exposure  in  the  same  experiment. 

The  impact  of  dry  gas  concentrations  of  1,  2,  and  10  ppm  CO  on 
the  cell  performance  over  a  wide  range  of  operating  conditions  was 
studied  by  combining  the  developed  methodology  with  a  previously 
described  high-resolution  gas  analysis  technique.  Over  the 
concentration  range  studied,  the  greatest  performance  impact  was 
found  to  occur  at  high  contaminant  concentration,  low  cell 
temperature,  low  anode  relative  humidity,  and  high  current  density. 
While  cell  temperature  had  the  most  effect,  the  self-induced 
performance  recovery  from  the  steady  state  poisoning  state  oper¬ 
ating  conditions  was  found  to  be  much  less  dependent  on  operating 
conditions  or  contaminant  concentration.  Lower  temperatures 
resulted  in  longer  recovery  times  although  the  final  performance 
recovery  was  similar  regardless  of  operating  conditions  and  led  to 
nearly  complete  recovery  after  approximately  5  h  of  recovery  time. 

Comparison  of  CO  exposure  experiments  with  different  cathode 
reactants  of  hydrogen,  air,  and  oxygen  provided  strong  evidence 
that  oxygen  permeating  through  the  membrane  from  the  cathode  to 
the  anode  significantly  reduced  the  impact  of  the  CO  on  cell 
performance.  The  crossover  oxygen  was  believed  to  enhance  the 
conversion  process  of  CO  to  CO2.  The  conversion  of  CO  to  CO2  at 
steady  state  decreased  with  increased  operating  temperature  of  the 
fuel  cell.  This  may  have  been  related  to  increased  oxygen  crossover 
due  to  increased  amounts  of  liquid  water  in  the  cell  at  lower 
temperatures. 

This  work  provided  evidence  that  the  complexity  of  the  CO 
contamination  process  increased  for  low  CO  concentration  levels 
with  respect  to  previous  studies  at  higher  concentrations.  Oxygen 
crossover,  the  presence  of  liquid  water,  the  thickness  of  the 
membrane,  the  catalyst  loading  and  transient  operating  conditions 
will  have  to  be  accounted  for  in  future  hydrogen  fuel  quality 
specifications.  Additional  model  development  and  experiments  at 
the  fuel  purity  recommendations  from  the  Society  of  Automotive 
Engineers  (SAE)  and  the  International  Organization  for  Standard¬ 
ization  (ISO)  working  committees  are  needed.  Future  work  at  our 
laboratory  will  address  the  experimental  aspects  of  these  needs 
while  Appendix  B  provides  numerical  fits  to  the  data  presented  in 
this  paper  which  we  hope  will  provide  modelers  access  to  valuable 
data  for  further  elucidation  of  CO  contamination  processes. 
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Appendix  A.  Methodology  for  analysis  of  cell  voltage 
response 

To  determine  the  range  of  cell  performance  degradation  for 
PEMFCs  operated  with  neat  H2/air  over  the  maximum  length  of 


a  contamination  experiment  several  MEAs  from  the  same  produc¬ 
tion  batch  were  tested  at  various  operating  conditions.  The  test 
duration  was  200  h,  which  was  the  maximum  experimental  time  of 
this  work.  The  results  of  three  experiments  are  shown  in  Fig.  10. 
Each  of  the  plots  represents  the  performance  of  a  different  MEA 
over  time  with  the  exception  of  the  plot  where  the  performance 
measurement  was  interrupted.  In  this  case,  the  cell  was  shut  down 
at  100  h,  then  restarted  and  operated  for  an  additional  100  h  to 
measure  the  performance  degradation  before  and  after  shutdown. 
The  operating  conditions  are  given  within  the  figures. 

Also  shown  in  Fig.  10  are  the  performance  degradation  rates 
calculated  from  a  linear  regression  of  the  data.  For  the  cell  shut¬ 
down  after  100  h  the  linear  degradation  rate  was  240  pV  h-1  before 
and  after  shutdown  indicating  that  the  cell  could  be  stopped  and 
restarted  without  affecting  the  degradation  rate.  The  linear 
degradation  rate  for  the  other  MEA  operated  at  80  °C  and  100  % 
relative  humidity  was  51  pVh-1.  This  indicated  that  a  range  of 
degradation  rates  could  be  observed  for  MEAs  from  the  same 
production  batch  operated  at  identical  conditions.  The  MEA  oper¬ 
ated  at  60° C  had  a  degradation  rate  of  125  pVh-1.  This  suggested 
that  the  MEA  degradation  rates  were  dependent  on  the  MEA  tested 
and  not  significantly  affected  by  the  cell  temperatures  used  in  this 
work.  Because  of  the  wide  range  of  observed  degradation  rates  it 
became  apparent  that  degradation  rates  needed  to  be  determined 
in-situ,  i.e.  during  the  contamination  experiment,  for  each  MEA  and 
operating  condition. 

Several  similar  200  h  experiments  were  analyzed  to  determine 
the  fraction  of  time  of  the  experiment  required  to  accurately 
extrapolate  cell  performance  over  the  entire  duration  of  the 
experiment.  The  data  were  fit  over  different  fractions  of  time  using 
Eq.  (3).  All  fits  started  at  the  beginning  of  the  experiment  and 
ranged  from  1  /6  to  the  entire  experiment  time,  e.g.  from  33  to  200  h 
for  a  200  h  experiment.  Fits  were  subsequently  extrapolated  over 
the  entire  duration  of  the  test. 

Fig.  11  shows  a  representative  example  experiment  performed 
at  80  °C  and  an  anode  inlet  humidity  of  80%.  The  figure  shows 
a  comparison  of  the  cell  voltage  data  with  data  extrapolated  from 
the  different  fitting  ranges.  The  voltage  data  fit  to  the  entire  dura¬ 
tion  of  the  experiment  showed  excellent  agreement  with  the 
measured  voltage  and  emphasized  the  accuracy  of  the  fitting 
equation.  As  the  amount  of  data  used  for  fitting,  i.e.  the  fitting 
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An/Cath:  80°C,  48.3/48.3  kPag,  2/2stoich,  80/50%RH,  H^Air 


0  50  100  150  200 

Time  [hours] 

Fig.  11.  Extrapolations  of  data  fit  with  Eq.  (3)  from  different  time  intervals  of  the 
experiment  for  cell  operated  for  200  h  at  80  °C  and  an  anode/cathode  inlet  relative 
humidity  of  80/50%.  All  fits  start  from  t  =  0. 


30  50  66.7  100  200 

Upper  Fitting  Limit  [hrs] 

Fig.  12.  Mean  square  error  associated  with  extrapolation  from  fits  of  the  five  time 
intervals  in  Fig.  13  for  cell  operated  for  200  h  at  80  °C  and  an  anode/cathode  inlet 
relative  humidity  of  80/50%. 


Table  4 

Operating  conditions,  contaminant  exposure  time  ranges,  fitting  time  ranges,  employed  fitting  equations,  and  fitting  parameter  and  R2  values  of  conducted  experiments  and 
fits. 


An/Ca 

conditions 

Current 
density 
[A  cm-2] 

CO 

concentration 

[ppm] 

Poisoning 
time 
range  [h] 

Fitting 

equation 

Fitting  range  [h] 

Contamination 

1 

Recovery 

Pi 

P2 

P3  P4 

P  5 

R2 

60  °C  48.3/ 

1 

1 

45.27-100.27 

17 

8.6-100.27 

-0.00025 

0.17308 

59.22434  17.37171 

0.99918 

48.3  kPag2/2 

22 

0-45 

0.00601 

0.05984 

-0.27753  -0.81608 

0a 

0.99852 

stoich 

1 

2 

40.27-90.27 

17 

4-90.27 

-0.0002 

0.2453 

46.91563  36.08006 

0.99892 

100/50%RH 

22 

0-40 

0.00333 

0.04822 

-0.11465  -0.74199 

0a 

0.9977 

1 

10 

40.4-70.49 

18 

39.27-41.2 

0.00005 

4.82588 

39.77731  0.56744 

0.99846 

19 

41.14-50 

0.27662 

-0.00435 

41.12917  0.04305 

0.99653 

20 

49-70.49 

0.28428 

9.77387E-4 

0.91612 

22 

0-40 

-0.0071 

0.08584 

-0.04685  -0.47243 

0a 

0.99318 

80  °C  48.3/ 

1 

1 

50.27-130.27 

17 

22.5-130.27 

-0.00017 

0.02566 

62.39081  14.27875 

0.98897 

48.3  kPag2/2 

22 

0-40 

0.0097 

0.0199 

-1.03109  -0.96593 

0a 

0.92548 

stoich 

1 

2 

40.27-90.27 

17 

6.5-90.27 

-0.00011 

0.04463 

45.76363  27.5232 

0.99855 

100/50%RH 

22 

0-26.1 

-0.00216 

0.0194 

-0.13905  -0.47695 

0a 

0.96912 

1 

10 

40.27-70.27 

18 

4-42.17 

0.00006 

5.80598 

39.68907  -1.05889 

0.99652 

19 

42.15-70.27 

0.20752 

-0.00007 

42.11886  0.01257 

0.9635 

22 

0-40 

0.00664 

0.01848 

-0.04652  -0.82592 

0a 

0.98442 

80  °C  48.3/ 

1 

1 

45.27-100.27 

17 

1-100.27 

-0.00012 

0.02473 

53.85678  18.83773 

0.99782 

48.3  l<Pag2/2 

22 

0-45 

-0.00078 

0.01825 

-0.65148  -0.63686 

0a 

0.95346 

stoich 

1 

2 

40.27-90.27 

17 

2-90.27 

-0.00013 

0.04853 

46.63424  24.76019 

0.99757 

80/50%RH 

22 

0-40 

0.00945 

0.0186 

-0.52559  -1.19436 

0a 

0.9555 

1 

10 

40.38-70.48 

18 

4-41.92 

0.00003 

7.34888 

39.98319  -1.60282 

0.99619 

19 

41.92-52 

0.18768 

-0.00059 

41.89154  0.02101 

0.96391 

20 

52-70.48 

0.19123 

7.52365E-4 

0.88771 

22 

0-40 

0.01023 

0.02192 

-0.02485  -0.65766 

0.00035 

0.99403 

80  °C  48.3/ 

1 

1 

50.27-130.27 

18 

27.3-56 

0.00032 

1.01504 

51.73667  -0.09594 

0.89508 

48.3  kPag2/2 

17 

53-130.27 

-0.01211 

0.02772 

59.30433  6.59218 

0.98808 

stoich 

22 

0-30 

0.00812 

0.1905 

-3.70526  -1.75432 

0a 

0.98852 

60/50%RH 

1 

2 

40.27-70.27 

17 

16-70.27 

-0.00028 

0.07737 

46.33266  25.75828 

0.99896 

22 

0-40 

0.00524 

0.02661 

-0.19723  -0.65537 

0a 

0.99405 

1 

10 

28.94-58.94 

18 

28.94-30.5 

0.00003 

7.11219 

28.39885  -1.44645 

0.99725 

19 

30.5-58.94 

0.26882 

-0.00049 

30.49458  0.02182 

0.98302 

22 

0-40 

-0.00085 

0.01974 

-0.03141  -0.80609 

0a 

0.99698 

60  °C  48.3/ 

0.2 

2 

85.18-215.18 

17 

5-215.18 

0.00006 

0.01876 

117.99301  11.69972 

0.99593 

48.3  l<Pag2/2 

0.6 

2 

50.38-120.6 

20 

5-51.65 

0 

0 

stoich 

21 

51.65-68.6 

0.02164 

0.00493 

58.48721  0.00016 

2.1309E-6 

0.99965 

100/50%RH 

17 

68.6-120.6 

-2 

0.15948 

53.90687  155,733 

0.93051 

22 

0-30 

0.00393 

0.05023 

-0.24847  -0.83531 

0.99907 

Forced  value. 
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Fig.  13.  Fitting  example  for  contaminant  experiment  with  dry  gas  concentration  of 
10  ppm  CO  at  80  °C  with  an  anode/cathode  inlet  relative  humidity  of  100/50%.  The 
example  data  was  partitioned  into  three  fitting  sections  covering:  (i)  the  pre-poisoning 
phase  and  the  beginning  of  transient  response  of  the  poisoning  phase,  (ii)  the 
remainder  of  the  poisoning  phase  including  the  steady  state  poisoning  state,  and  (iii) 
the  recovery  phase. 


y  =  p-i+P2-x  (is) 

y  =  Pl  +P2-(X-P3)+P4-(X-P3)2+P5-(X-P3)4  (16) 

y  =  Pl  +  (P2'(X-p3)P4)+P5-X  (17) 

Table  4  lists  the  operating  conditions,  contaminant  exposure 
time  ranges,  fitting  time  ranges,  employed  fitting  equations,  and 
fitting  parameter  values  of  all  conducted  experiments,  and  the  R 2 
value  of  the  fits.  With  this  information  the  A r/  versus  time  of  each 
experiment  can  be  accurately  reconstructed.  Note  that  (i)  data 
values  prior  to  contaminant  exposure  can  be  set  to  zero,  and  (ii) 
extrapolation  of  the  fits  beyond  the  indicated  fitting  ranges  will 
likely  result  in  significant  deviation  from  the  measured  values. 
Fig.  13  shows  an  example  for  fitting  the  A rj  of  an  experiment  using 
three  equations:  (i)  the  exponential  Eq.  (13)  was  employed  to 
describe  the  pre-poisoning  phase  and  the  beginning  of  the  non¬ 
steady  state  ranging  from  4  to  42.17  h,  (ii)  Eq.  (14)  which  con¬ 
tained  linear  and  logarithmic  terms  was  employed  to  describe  the 
remainder  of  the  non-steady  state  and  the  steady  poisoning  state 
period  ranging  from  42.17  to  70.27  h,  and  (iii)  Eq.  (17)  which  con¬ 
tained  linear  and  power  terms  was  used  for  cell  performance 
recovery  ranging  from  70.27  to  100.27  h.  The  recovery  time  is 
shown  as  a  second  time  scale  in  Fig.  13  which  starts  at  a  time  of 
zero.  This  time  scale  should  be  used  when  reproducing  the  data  and 
was  employed  to  maximize  the  fitting  accuracy  of  the  cell  perfor¬ 
mance  recovery. 


range,  was  decreased,  however,  the  accuracy  of  the  extrapolation 
decreased. 

To  determine  the  fraction  of  experiment  time  necessary  to 
provide  a  reasonably  accurate  extrapolation  the  mean  square  error 
(MSE)  of  the  measured  voltage  data  to  the  data  extrapolated  from 
the  fit  was  calculated  for  each  fitting  range  using  Eq.  (11).  Fig.  12 
shows  the  MSEs  for  the  data  and  fits  shown  in  Fig.  11.  The  data 
indicated  that  the  MSE  of  the  extrapolated  data  significantly 
improved  when  the  fitting  range  was  at  least  one  third  of  the  entire 
test  duration. 

M£E  ELi[Wtj)-VFit(ti)]2  (11) 

n  1  ; 

Appendix  B.  Mathematical  data  description  for  modeling 
efforts 

Research  institutions  throughout  the  world  are  working  on  the 
impact  of  contaminants  on  fuel  cell  performance  and  many 
contribute  to  national  and  international  efforts  to  develop  a  purity 
standard  for  hydrogen  used  in  PEMFC  applications,  specifically 
PEMFC  vehicle  applications.  Low-level  contaminant  experiments 
and  contaminant  model  development  are  essential  for  the  success 
of  this  effort.  In  this  section  we  report  our  data  to  the  research 
community  to  make  it  available  for  model  development  efforts. 
The  overpotential  change  of  the  cell,  Ary,  for  each  experiment  is 
reported  using  the  mathematical  expressions  Eqs.  (12)-(17). 
While  these  equations  and  their  fitting  parameters  do  not  contain 
any  physical  meaning,  they  served  as  a  convenient  medium  to 
transfer  the  large  amount  of  experimental  data  presented  and 
discussed  in  this  work. 


y  =  P2  +  (Pl  -P2)/(l  +  (*/P3)P4) 

(12) 

y  =  Pi-exp(p2-(x-p3)+p4-(x-p3)2) 

(13) 

y  =  Pi  +P2-(x-P3)+P4-ln(x-p3) 

(14) 

Glossary 

Symbols 

MSE  mean  square  error 

n  number  of  data  points 

Pi,  P2,  P3,  P4  fitting  parameter 
QcBai  quality  of  carbon  balance 

tarb  arbitrary  chosen  saturation  time 

t  time 

tSat  time  required  to  reach  saturation  of  catalyst  layer  with 
contaminant 

tss  time  steady  state  poisoning  state  is  reached 

tstop  time  contaminant  exposure  stops 

1/ceii  cell  voltage 

Vfit  fitted  voltage 

At/  overpotential  change 

?co  magnitude  of  CO  conversion 
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